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Modified  Linear  Theory  F-function  Applied  to  Numerical 

Prediction  of*  Sonic  Booms 

Timothy  R  Jung  1  ,  Ryan  R  Starkey  2  and  Brian  Argrow  3 
University  of  Colorado ,  Boulder,  CO ,  80309,  USA 

Sonic  boom  pressure  profiles  calculated  by  modified  linear  theory  are  evaluated-  These  methods 
transform  parametric  description  of  aircraft  into  F-functions,  which  are  used  to  predict  the  pressure 
profile  after  propagation  through  the  atmosphere*  Although  Whitham’s  F-function  equation  is  the 
most  widely  used  method,  it  cannot  be  applied  to  aircraft  that  have  discontinuities  in  their  area  pro¬ 
files-  Thus,  Light  hill's  F-function  is  shown  to  be  the  best  method*  Procedures  are  presented  to  adapt 
LighthilPs  method  so  it  can  he  used  for  numerical  calculation*  Then,  sonic  boom  pressure  profiles  from 
wind  tunnel  experiments,  computational  fluid  dynamics,  and  flight  tests  are  compared  to  profiles  from 
Lighting  method  to  explore  its  accuracy*  In  general,  LighthilPs  modified  linear  theory  method  is  able 
to  predict  the  sonic  boom  shape,  magnitude  and  duration  within  U\%  and  account  for  variations  in 
lift,  Mach  number  and  propagation  angle-  The  final  case  investigated  is  the  Its  complex  shape 

is  not  easily  transformed  into  an  axi-symmetric  area  profile,  reducing  the  accuracy  of  the  modified 
linear  theory  pressure  profile-  However,  this  study  demonstrates  that  modified  linear  theory  is  a  fast 
and  sufficiently  accurate  to  predict  sonic  booms  and  is  appropriate  for  the  preliminary  and  conceptual 
design  stages  of  quiet  sonic  boom  aircraft* 


Nomenclature 

co  =  drag  coefficient 

cl  =  lift  coefficient 
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/  =  source  strength 

F  -  F-function,  \/W 

h  =  function  in  LighthilFs  F-function 

=  propagation  distance  in  vertical  direction,  ft 
H  -  Heaviside  step  function 

ISPR  —  initial  shock  pressure  rise,  Jb/ft2 
K±  =  Bessel  function  of  second  kind 

C  =  length  of  aircraft,  ft 

L  =  lift,  lbs 

M  =  Mach  number 

p  ~  Heavi  s  ide 7s  opera  tor 

poo  ~  freestream  static  pressure,  Ib/ft2 

Ap  =  pressure  perturbation,  Ib/ft2 

Ap  =  non-dimensional  pressure 

£  -  dynamic  pressure,  Ib/ft2 

r  -  radial  coordinate,  ft 

R  -  radius  of  body  of  revolution,  ft 

Sy  =  cross-sectional  area,  ft2 

Sl  =  cumulative  lift  along  x  axis  in  units  of  area,  ft2 

Seq  =  total  equivalent  area,  ft2 

u  “  velocity  perturbation  in  x  direction,  ft/sec 

v  =  velocity  perturbation  in  r  direction,  ft/sec 

W  =  weighting  factor  in  F-function  ft-1/2 

x  =  longitudinal  coordinate  from  nose  of  aircraft,  ft 

x  =  non-d i  m  en si  o  nal  x  d  i  stan  ce 

Ax  =  step  size  in  x,  y  direction,  ft 

y  =  ch  aracteri  st  i  c  coord  i  n  ate,  ft 

fi  =  Prandtl-Glauert  coefficient,  \JM2  -  1 
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7  =  rati  o  o  f  spec  l  fi  c  h  eat  s 

C  ~  du  mmy  vari  ah le  of  i  n  tegrat  ion  along  y 

y.  =  Mach  angle,  radians 

X  =  dummy  argument 

4>  =  velocity  potential 

=  propagation  angle,  degrees 


Subscripts 

eq 

=  total  equivalent  area 

f 

=  fuselage 

n 

=  summation  index 

i 

—  location  i 

L 

=  Ligbthill 

ML 

=  modified  Lighthill 

V 

=  volume 

W 

=  Whltham 

L  Introduetkm 

The  University  of  Colorado  has  used  conceptual  and  preliminary  design  software  named  RapidP  to  de¬ 
sign  quiet  sonic  boom  aircraft  [  1  ]  and  to  investigate  how  to  conduct  scaled  sonic  boom  tests  in  the  atmosphere 
using  Unmanned  Aerial  Systems  [2].  Since  RapidF  is  used  as  a  conceptual  design  tool,  it  must  be  fast,  and 
thus,  lower  fidelity;  however,  it  must  also  be  accurate  so  that  the  designer  can  property  evaluate  a  config¬ 
uration,  Therefore,  modified  linear  theory  was  selected  as  the  foundation  of  RaptdF  because  it  is  capable 
of  quickly  and  accurately  predict  far  field  pressure  profiles,  A  pictorial  representation  of  RapidF’s  flow  is 
shown  in  Fig.  J .  The  program  accepts  design  inputs  and  flight  parameters.  The  most  important  inputs  are  are 
shown  in  Table  1  These  values  are  transformed  into  an  equivalent  a  xi-sym  metric  area  and  lift  profile  along 
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ihe  aircraft  s  longitudinal  axis,  as  shown  in  the  top  of  the  Rapid!7  box  in  Fig.  I.  Then,  using  these  area  pro¬ 
files,  RapidF  calculates  the  F-function  using  modified  linear  theory.  Finally,  RapidF  sends  the  F-function  to 
CUBoom  w  hich  is  a  geometric  acoustic  propagation  subroutine  to  predict  the  design's  sonic  boom  pressure 
profile  on  the  ground. 

For  brevity,  the  scope  of  this  paper  is  limited  to  the  F-function  calculations  and  resultant  pressure  profiles. 
Thus,  discussion  concerning  the  methods  used  to  generate  the  area  profiles  is  abridged,  A  comprehensive 
description  of  RapidF's  methodology  can  be  found  in  Ref.  [3f  The  objectives  of  this  paper  are:  1)  to 
determine  the  best  F-function  equation  2)  specify  how  to  use  this  equation,  and  3)  demonstrate  that  modified 
linear  theory  is  as  nearly  as  accurate  as  CFD  in  predicting  a  design’s  ground  pressure  profile,  while  taking 
significantly  less  time  to  execute. 


Fig.  1  Overview  of  RapidF 


II.  Modified  Linear  Theory 

Modified  linear  theory  is  based  on  axi-symmetric  linear  theory,  as  the  name  implies.  By  limiting  the 
flow  to  being  above  transonic  but  below  hypersonic  speed;  and  assuming  a  steady,  irrotational  and  isentropic 


flow  field  with  small  perturbations,  linear  theory  reduces  the  continuity,  momentum  and  energy  equations  to 
a  single  equation  with  one  variable,  the  velocity  potential  [4,  5]: 


d2<t>  1  M  _  ,2^  _  n 

dr2  r  dr  dx2 


(1) 
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Table  I  Summary  of  inputs  into  RapidF  software 


Fuselage 

Surfaces 

Engine 

Flight  Conditions 

Diameter 

Location 

Location 

Mach  number 

Length 

Leading  edge  sweeps  (2) 

Length 

Angle  of  Attack 

Diameter  changes  Trailing  edge  sweeps  (2) 

Radius 

Altitude 

Nose  type 

Chord  and  span 

Internal /external 

Weight 

Nose  length 

Airfoil  shape 

Pressure  ratio 

Boat  tad  type 

Incidence  angle 

%  Spillage 

Boat  tail  length 

Dihedral 

Fuel  Consumption 

Twist 

Thrust 

where  the  PnmdtLGlauert  coefficient  is: 


0  =  = 


1 

tan  fi 


(2) 


As  shown  in  Fig.  2,  the  characteristics  are  linear  of  the  form  y  =  x  -  fir.  Linear  theory  has  been  widely 
used  to  obtain  first-order  approximations  of  surface  pressures  and  aerodynamic  forces  on  aircraft  and  has 
consistently  compared  well  to  experimental  data  [6],  However,  as  presented  above,  it  cannot  predict  the  flow 
field  lar  from  a  supersonic  source,  which  is  the  location  of  concern  for  sonic  boom  minimization  research, 
because  the  theory  fails  to  account  for  non-linear  effects. 

Modified  linear  Lheory,  pioneered  by  Whitham  [7],  adapted  linear  theory  to  be  used  away  from  the 
source.  He  used  a  potential  defined  by  a  source  function,  /(at),  along  the  a:- ax  is  as  a  solution  to  Eq.  (1): 

,.T-13r  -f(c) 

9~  .L  v^'-O2  -/32r2<iC  (3) 


The  source  strength  is  determined  by  the  boundary  conditions.  As  with  any  potential  function,  the  velocity 
profile  can  be  found  by  taking  the  derivative  of  Eq.  (3)  with  respect  to  x  and  r.  Whitham  used  an  F-function: 

c  ,4, 

Jo  vv-  C 

to  simplify  calculations  because  the  F-funetion  remains  constant  along  a  characteristic  and  the  flow  properties 
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Fig*  2  Linear  theory  geometry. 


can  be  easily  calculated  via: 


dtp 

dx 


(5) 


V  £T 


A p  = 


AT2 


(7) 


In  Eq.  (4),  the  integration  is  conducted  over  y.  As  shown  in  Fig.  2,  pi  can  be  interpreted  as  the  location  on 
the  .x-axis  where  the  characteristic  caused  by  the  radius  change  at  x\  starts.  Thus,  integration  over  y  is  in  the 
-direction,  along  the  aircraft  centerline*  Once  the  F-function  is  calculated  from  Eq.  (4),  the  complete  flow 
field  can  be  estimated  Irom  Eqs,  (5-7).  The  velocities  are  only  depended  on  the  characteristic,  y  and  radial 
distance,  r.  In  general,  as  can  be  seen  by  inspecting  Eqs.  (5-7)  the  perturbations  decrease  as  a  function  of 
r~1/2.  This  also  causes  the  perturbations  to  be  infinite  at  the  source,  limiting  modified  linear  theoiy  to  the 
region  away  from  the  source  or  near  the  nose  where  fir  3>  y. 

Equations  (5—7)  are  not  sufficient  in  ol  themselves  to  describe  the  far  held,  because  variations  in  the 
local  speed  of  sound  will  cause  the  characteristics  to  curve.  Whitham  modified  the  characteristics  to  include 
a  non-linear  term  to  allow  them  to  bend  [7]: 


(8) 
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The  modified  charge  ten  sties  can  be  used  to  predict  shock  locations.  As  shown  in  Fig.  3,  each  characteristic 
curves  at  its  own  rate  related  to  its  F- function  value,  which  allows  characteristics  to  intersect.  When  two 
characteristics  merge,  a  shock  forms.  This  is  an  improvement  from  classical  linear  theory  which  cannot 
predict  shock  locations,  even  on  the  object's  surface.  By  assuming  the  shocks  are  weak  (small  entropy 
change),  the  shock  strength  can  be  calculated  from  Eq .  (7)  in  the  form: 

In  summary,  modified  linear  theory,  as  a  subset  of  linear  theory,  has  the  same  assumptions  and  restric- 
dons.  Its  prediction  of  the  flow  properties  is  still  first  order,  with  the  exception  of  the  characteristics,  which 
are  higher  order.  This  modest  change  enables  modified  linear  theory  to  be  used  in  the  far  field.  The  F- 
funciion  is  the  key  to  modified  linear  theory  because  it  has  multiple  applications.  In  Eqs.  (5-7),  it  is  used 
to  describe  the  velocity  and  pressure  perturbations  and  in  Eq.  (8),  it  describes  the  curvature  of  the  charac¬ 
teristics,  which  predicts  shock  locations  and  strengths.  For  this  reason  modified  linear  theory  is  well  suited 
for  RapidF.  When  designing  a  quiet  supersonic  aircraft,  the  user  is  immediately  aware  of  the  component's 
impact  on  the  F- function  and  thus,  the  sonic  boom  pressure  profile  as  was  demonstrated  in  Ref.  m. 


Fig,  3  Example  of  modified  linear  theory  curved  characteristics  from  an  elliptical  body. 
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A*  Whitham’s  F-fmieticm 


By  using  the  How  Uingeney  boundary  condition,  Whitham  determined  a  method  to  calculate  the  F- 
function  from  the  surface  shape.  The  radial  velocity  can  be  found  by  taking  the  derivative  of  Ecj.  (3)  with 
respect  to  r: 


-«-'r  (» - o/'(o 

dr  r  Jo  y/{x  -  02  -  fpr2 

Evaluating  the  above  equation  at  the  surface  and  assuming  a  slender  body  (R  s=  0): 


(10) 


|  body  — 


f(y) 


(ID 


Since  the  (low  must  be  tangeni  to  the  body  at  the  surface: 


R'  = 


r.1  \hody 

1  ""f"  tr  1 b,  Ilf 


‘  |  body 


(12) 


Combining  the  above  two  equations,  the  source  strength  can  be  calculated  from  the  surface  geometry  by: 


f(y)  =  RR' 

The  cross  sectional  area,  Sy<  associated  with  R  is: 


(13) 


and  by  differentiation: 


Sv  =  7T  R“ 


(14) 


Sy  =2nRR' 


(15) 


Combining  Eqs,  ( 1 3)  and  ( 1 5): 


f(y) 


(16) 


Finally,  by  differentiation: 


f(y)  =  ^s'v(z)  (1?) 

which  can  be  substituted  into  Eq.  (4)  to  provides  a  method  for  calculating  the  F-function: 
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Equation  (18)  will  be  referred  to  hereafter  as  Whitham’s  F- function  equation.  Results  from  this  equation  can 
be  used  in  Fqs.  (5-7)  to  describe  the  How  field.  However,  it  requires  Sfv  and  Sy  to  be  continuous. 

Whitham  initially  limited  his  theory  to  axi-sym  metric,  non-lifting  objects.  Later,  he  [8]  and  Walkden 
[9]  expanded  the  theory  to  include  non -axi -symmetric  and  lifting  surfaces,  respectfully,  which  were  experi¬ 
mentally  validated  by  Carlson  [10],  In  general,  these  researchers  added  the  equivalent  area  due  to  lift  to  the 
equivalent  area  from  volume  to  produce  the  total  equivalent  area: 

Scq{y)  =  SL(y)  +  Sv(y)  (19) 

where  the  area  due  to  lift,  is  the  cumulative  lift  in  units  of  area: 

Si(y)  =  — ^  L{Qd(  (20) 

The  total  equivalent  area  (referred  to  as  total  area  for  brevity),  Seq(y),  is  the  sum  of  the  two  areas  and  is  used 
in  Eq,  ( 1 8)  to  calculate  the  F-function.  Thus,  Whitham *s  F- function  can  be  used  on  a  no n- ax i -symmetric 
body  producing  lift. 

To  find  the  F-function,  using  Whitham’s  F-function  at  location  y,  the  second  derivative  at  all  locations  Q 
in  front  of  y  are  weighted  by: 

*W0  =  y=  (21) 

Based  on  this  weighting,  one  can  interpret  Eq.  (18)  as  stating  (hat  the  F-function  is  mostly  a  function  of 
the  local  curvature.  The  curvature  ol  locations  in  front  of  the  point  of  interest  do  influence  the  F-function; 
however,  the  influence  rapidly  decreases  as  the  distance  between  the  points  increases.  For  example,  Fig.  4(b) 
compares  the  local  first  and  second  derivative  to  Che  F-function  of  a  low  sonic  boom  aircraft  from  Ref. 
[  I  ].  Since,  Whitham’s  F-function  integrates  the  second  derivative,  the  first  derivative  follows  the  F-functions 
trends;  thus,  the  first  derivative  can  be  used  as  a  crude  approximation  of  the  F-function.  This  is  consistent  with 
supersonic  linear  theory  where  the  surface  pressure  is  only  a  function  of  the  local  surface  inclination  and  that 
information  cannot  be  transmitted  upstream.  However,  Fig.  4  highlight’s  a  problem  when  using  Whitham’s 
F-function  on  an  aircraft.  Whitham’s  F-function  equation  requires  S"  to  be  continuous.  Although  aircraft 
are  usually  smooth  to  minimize  drag,  discontinuities  do  exist  in  the  cross-sectional  area  profile,  such  as  the 
engine  inlets,  as  shown  in  Fig.  4(a).  Whitham’s  F-function  in  Fig.  4(c)  shows  a  large  spike  in  the  F-function 
at  this  location.  These  spikes  are  sensitive  to  grid  density  as  will  be  shown  later. 
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(a)  Area  profiles. 


(b)  F- function  and  derivatives  of  area. 


(e)  Aircraf  profile  view  and  F-function, 

Fig.  4  Area  profiles,  derivatives,  and  Whit  ham’s  F-funetlon  of  an  aircraft* 


To  use  Whitham’s  equation,  in  RapidF,  it  must  be  discretized  to  allow  numerical  integration.  Using  a 
second  order  estimate  of  the  second  derivative  of  the  area,  Eq,  (4)  becomes: 


Pi 


~  2 Sn  -r  $n  +  l  3. 

{ Ax ) 


Ax 


(22) 


Integration  cannot  include  n  —  i  because  the  last  term  would  be  infinite.  For  the  last  term  in  the  summation. 


n  =  ?'-  !;  thus,  x*  -  xn  —  Ax  and  the  term  simplifies  to: 


^  (Si^2  ~  2Sj_i  +  Si)  Ax~3/2  (23) 

Thus,  as  the  user  increases  the  grid  density,  decreasing  Ax,  the  weight  on  this  last  point  increases,  making 
Whitham's  equation  sensitive  to  grid  density.  Because  of  this  issue,  combined  with  the  inability  to  account 
for  discontinuities,  a  more  robust  equation  was  sought  for  RapidF. 


10 


B,  Stieltjes  F- function 


One  method  to  calculate  an  F-function  with  discontinuities  is  to  use  a  Stieltjes  integral  of  the  general 
form  [1 !]; 


I  P(x)dG(x )  =  J  P{x)G'{x)dx 

Thus,  Whitham’s  F-function  can  be  calculated  by  [5]: 


(24) 


y  dS\ 0 
Vy  -  C 


(25) 


This  equation  will  be  referred  to  as  the  Stieltjes’  F-function,  It  still  requires  a  correction  for  where  S,(C)  is 


discontinuous.  According  to  Whitham  [7],  this  correction  is: 


Am‘  =  i]fw,As:  m 

which  is  essentially  the  same  as  Lighthill’s  F-function  (discussed  in  the  next  section).  The  above  equations 
can  be  applied  in  the  following  manner:  Equation  (25)  is  used  to  calculate  the  F-function  up  to  a  discontinuity 
in  S'(C),  then  Eq.  (26)  is  used  to  calculate  the  F-function  jump  at  the  discontinuity.  Then,  Eq.  (25)  is  used 
again  until  the  next  discontinuity.  The  process  is  repeated  until  reaching  the  end  of  the  aircraft.  However, 
Eq.  (25)  must  be  calculated  numerically.  To  do  so,  dS '  at  point  n  is: 


dSL 


C f  _  Q t 

^ forward  ^backwards 

Sn+ 1  ~  Sn  ^  Sn  —  Sn  ^  j 

Ax  Ax 

Sn-  \  ~  2Sn  H~  ffre+1 

Ax 


Substituting  this  into  Eq.  (25)  and  discretizing: 


Tl—l  '■ 


_  i  —  2Sn  +  Sn  + 1 


(27) 


Ax 

This  is  numerically  the  same  as  Eq.  (22)  and  there  is  no  advantage  of  the  Stieltjes  integral  compared  to 
Whitham’s  F-function.  For  example,  an  F-function  calculated  from  Eq.  (25)  is  compared  to  one  from 
Whitham’s  Function,  Eq.  (18),  in  Fig.  5.  As  expected,  they  are  identical.  The  Stieltjes  F-function  had  to 


be  offset  vertical  because  the  F-functions  are  virtually  the  same.  This  observation  has  been  previously  stated 
by  Ritzel  and  Gottlieb  in  Ref.  [5],  Therefore,  the  Stieltjes  F-function  equation  was  not  pursued  as  a  viable 
choice  for  RapidF. 


II 


Fig,  5  Whit  ham's  and  Stieltjes*  F- functions  (Stieltjes  offset  by  0.1  ftI/2), 

C,  LighthilFs  F-function 

In  his  1952  paper  [7],  Whithatn  built  on  Lighth ill’s  linear  theory  for  bodies  with  slope  discontinuities 
[  12]  to  develop  an  alternative  to  Eq.  (18): 

F{,)=r1pr  ,/i^o^  {28) 

Jo  \j  mo  Amo )  ^  (  8) 

The  above  equation  will  be  referred  to  as  LighthilFs  F-function,  where: 

hM  =  '») 

Here  p  is  the  Heaviside  operator  for  differentiation,  K1  is  the  Bessel  function  of  second  kind,  and  H(x)  is  a 
step  (unction.  Hie  derivation  of  Lighth  ill's  E- (unction  did  not  require  the  assumption  of  a  smooth  area  profile 
as  was  in  Whithatrfs  and  Stieltjes"  P-functions,  This  is  a  significant  improvement,  which  allows  Lighth  ill's 
F-function  to  be  applied  to  objects  that  have  discontinuities.  For  example,  a  comparison  of  Whitham  and 
Ugh  thill  methods  is  shown  in  Fig.  6  for  a  nose  cone  to  fuselage  joint  with  a  3,2°  change  in  flow  angle. 
Although  the  discontinuity  is  small,  Whitham  predicts  a  large  expansion.  This  example  demonstrates  that 
even  a  small  discontinuity  has  significant  effect  on  WhithanFs  F-function.  Since  a  digitized  representation 
of  an  aircraft's  area  profile  has  a  finite  jump  at  each  grid  point  as  shown  in  Fig.  7,  any  numerical  area  profile 
effectively  has  a  discontinuity  at  every  location  and  LighthilFs  method  appears  better  suited  for  calculating 
an  aircraft  F-function  using  numerical  integration. 

The  only  difference  between  LighthilFs  and  WhithanFs  F-functions  is  the  weighting.  Since  LighthilFs 
equation,  Eq.  (28),  is  also  a  Stieltjes  integral,  his  uses  of  dS'(C)  is  equivalent  to  Whitham’s  S"(£)dC ,  There- 
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fore,  the  large  disparity  in  the  F-functions  shown  in  Fig.  6  is  not  caused  by  Whitham’s  use  of  6T"(C),  but  by 
his  weighting.  Based  on  Eq.  (28),  Lightill’s  F-function  weighs  each  point  by: 


Wl(  0 


(30) 


The  h(x)  function,  Eq.  (29),  is  plotted  in  Fig.  8.  When  *  <  -1.  then  h(\)  =  0;  when  x  =  -1,  then 
Mx)  —  1;  an£l  when  \  is  tefge.  h(x)  ~  (2x)_l^'%  and  the  weight  simplifies  to: 


WL 


mo 


mo  V  2  (v  -  0 


(31) 


and  comparing  this  to  Eq.  (21),  Whitham’s  and  LighthilJ’s  weighting  are  the  same.  Their  weights  are  only 


significantly  different  where  ( y  —  £)/(/3 /?(£))  is  small,  near  y,  as  shown  in  Fig.  9.  This  figure  demonstrates 


Fig.  6  Comparison  of  Whitham  and  Lighthill  F-functions  at  a  discontinuity.  [5]. 


(a)  Smooth  or  continuous.  (b)  Discrete. 

Fig.  7  Smooth  and  discrete  total  equivalent  area  profiles. 
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Fig*  8  The  h(x)  function* 


Fig,  9  Whilliam  and  Lighthill  weighting, 

that  at  large  distances  WhithanTs  and  LighthilTs  weighting  are  the  same  as  suggested  by  Eq.  (31),  It  is  also 
important  to  note  that  the  area  under  each  curve  is  equal. 

Light  hill  s  i  -function  has  not  been  widely  implemented  due  to  the  difficulty  in  implementing  Heavi¬ 
side  calculus  and  Bessel  functions.  In  the  1980s,  Gottlieb  and  Ritzel  compared  WhithanTs  and  LighthilFs 
methods  to  projectiles  with  discontinuities  [5,  1 3 1 *  They  found  that  WhithanTs  theory  poorly  predicted  the 
F- function  at  discontinuities,  which  should  have  been  anticipated  since  a  discontinuity  violates  the  smooth 
assumption.  One  drawback  of  using  Ligh thill's  method  is  a  higher  computational  cost.  In  order  to  speed 
processing  time,  Gottlieb  and  Ritzel  proposed  using  a  modified  method  where  WhithanTs  equation  is  used 
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in  regions  with  a  smooth  area  profile  and  adding  a  correction  at  discontinuities  based  on  LighthilVs  equation. 
Initially,  this  method  was  attempted  in  this  study,  but  was  difficult  to  implement  since  it  requires  logic  to  find 
discontinuities.  Furthermore,  every  location  is  essentially  a  discontinuity  when  doing  numerical  integration, 
as  was  shown  in  big.  7(b).  Although  Lighthilbs  F-function  is  valid  at  every  point,  Gottlieb  and  Ritzel  were 
hesitant  to  solely  use  LighthilTs  equation  because  of  computation  time  and  advocated  for  the  Stieltjes  method 
discussed  earlier,  in  the  20+  years  since  their  papers  were  published,  advances  in  computers  have  made  this 
argument  irrelevant.  Thus,  for  this  study  a  method  was  developed  to  implement  Eq,  (28)  at  every  location  for 
aircraft  designs* 

LighthilTs  method  is  more  computationally  expensive  than  WhithanTs  because  Eq.  (29)  is  time  con¬ 
suming  to  calculate*  Previously,  this  cost  was  avoided  by  approximating  h{x)  as  (2x)~1/2.  However,  if 
this  approximation  is  used,  Eq.  (28)  simplifies  to  Eq.  (18)  and  any  advantage  in  using  LighthilTs  method  is 
lost.  Thus,  in  order  to  improve  accuracy  and  computational  speed  a  look-up  table  was  created  for  k(x)>  as 
suggested  by  Gottlieb  and  Ritzel  [5].  This  allows  LighthilTs  method  to  be  used  accurately  and  with  no  speed 
penalty, 

LighthilTs  equation  has  never  been  applied  to  an  aircraft,  only  to  a xi- symmetric,  non-lifting  objects; 
therefore,  some  adaptations  needed  to  be  made  for  use  in  RapidF.  First,  the  total  equivalent  area  (lift  and 
volume)  from  Eq,  (19),  is  used  to  calculate  dSf{Q)  based  on  Walk  den  work  that  concluded  lift  could  be 
replaced  by  an  equivalent  area  defined  in  Eq,  (20)  [9]. 

Next,  the  domain  of  in  finance  needed  to  be  modified.  In  WhithanTs  F-function,  the  domain  of  influence 
is  limited  by  the  limits  of  integration.  However,  since  the  limits  of  Lighthill’s  equation  arc  zero  to  infinity, 
the  domain  of  influence  is  controlled  in  the  //  function  in  Eq.  (30),  The  ftR{(ft  in  the  denominator  of  the  input 
to  the  h  lunction  has  significant  effect  on  the  calculations,  which  is  shown  pictorially  in  Fig,  10.  As  seen  in 
Eq.  (30),  the  ft  R(  f)  term  provides  a  method  to  non-dimensionalize  y  -  £  and  also  determines  at  what  y  -£ 
value  the  h(X)  is  zero,  due  to  the  step  function.  In  other  words,  the  ftR{Q  controls  how  far  beyond  point  ?/, 
the  aircraft’s  shape  contributes  to  the  F-function  at  point  y,  as  shown  in  Fig.  10(b),  If  R(Q  is  approximated 
as  R(y).  then  the  domain  ol  influence  is  as  shown  in  Fig,  10(c),  which  is  consistent  with  linear  theory:  a 
radius  change  at  y  +  ft R  act  as  if  it  began  at  the  centerline  at  y. 

RapidF  slices  the  aircraft  at  the  Mach  angle  and  projects  the  area  along  characteristics  back  to  centerline 
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(a)  Points  in  front  of  y.  (b)  Points  after  y . 


tc)  Domain  of  influence  assuming  f?(<)  s;  R(y ). 

Fig*  JO  Lighthill’s  weighting* 

as  was  shown  in  Fig.  2.  For  example,  refer  to  the  engine  shown  in  Fig.  1  1(a).  It  is  located  at  127  ft  from  the 
aircraft  nose  and  is  7  ft  below  the  aircraft.  If  the  aircraft  is  flying  at  Mach  1.7,  then  the  area  increase  from 
the  engine  acts  as  if  it  emanates  from  1 17.6  ft  along  the  centerline.  Thus,  the  area  profile  from  the  engine 
is  shown  in  Fig.  11(b).  In  Hg.  1 1(c),  the  Light  hill  F-fu.nct  ion  from  this  area  profile  is  shown  compared  to 
Whitham  s  F- function.  Because  LighthiM  s  domain  of  influence  extends  past  the  point  of  interest,  as  was 
shown  in  Fig.  10(c),  the  engine  effects  in  the  LighthilFs  F-function  begin  at  approximately  111  ft,  where 
(111  11 7.6) /(/?#(£))  —  1.  This  is  an  an  undesirable  attribute  of  LiglithilFs  method  because  the  engine 
should  influence  the  F-function  at  1  17.6  It.  In  RapidF,  since  the  aircraft  has  already  been  sliced  at  the 
Mach  angle,  as  shown  in  Fig.  2,  the  equivalent  area  at  point  y\  includes  the  effects  of  the  area  and  lift  at 
X}  4-  0R,  Thus,  when  RapidF  calculates  the  total  equivalent  area,  it  has  already  projected  it  back  to  the 
aircraft  centerline,  consistent  with  the  domain  of  influence  shown  in  Fig,  10(c).  To  avoid  shifting  the  engine 
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{ a )  Geornct  r>r  of  t  he  engi  ne , 


fbj  Engine  area  profile. 


(c)  F-functions, 


(d)  Weights. 


Fig,  II  Engine  example. 


area  forward  a  second  time,  Lighlhiil’s  weighting  is  modified; 


(32) 


as  shown  in  Fig.  1 1(d),  with  the  corresponding  F-functimi  shown  in  Fig.  1 1(c).  This  allows  features  in  the 
F- function  to  appear  at  the  appropriate  location  in  the  F-function,  while  keeping  the  area  under  the  weighting 


curves  in  Fig.  1 1(d)  constant. 

Finally,  a  value  was  need  to  be  used  for  R(Q-  As  already  explained,  the  ft  R(Q  term  affects  the  effective 
domain  of  influence.  Several  values  were  explored  to  include: 


•  ^(0  =  \fsljir-.  Total  equivalent  area  (lift  and  volume),  in  units  of  length 
8  RviQ  —  s/SvJn:  Equivalent  area  of  volume,  in  units  of  length 


•  R/{ C):  Local  fuselage  radius 
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The  first  two  radii  choices  were  significantly  larger  than  the  fuselage  radius  and  shifted  the  effect  of  aircraft 


components  too  early  in  the  F-function.  Thus,  the  radius  used  in  RapidF  is  the  fuselage  radius.  However,  it  is 
possible  to  have  components  such  as  highly  swept  wings  and  horizontal  tails  to  extend  past  the  fuselage.  In 
these  cases  R/(()  =  0.  causing  a  divide  by  zero  problem  in  Eq.  (30).  This  is  avoided  by  setting  the  minimum 
Rf(0  vtdtie  to  0.02Rj<max.  All  of  the  above  modifications  to  Lighthill’s  F-function  result  in; 


This  equation  will  be  referred  to  as  the  modified  Lighthill’s  F-function  equation. 

D.  F-function  sensitivity 

A  grid  sensitivity  and  convergence  study  was  conducted  on  the  classical  Whitham  method,  Eq.  (18), 
and  the  modified  Lighihitl  method,  Eq.  (33).  The  primary  input  to  the  resolution  of  the  RapidF  solution  is 
the  user  defined  input  distance  between  the  slices  along  the  aircraft  centerline.  Ax.  However,  the  discussion 
below  presents  results  using  the  “Number  of  Points’"  or  grid  size,  which  is  the  aircraft  length  divided  by  Aar. 

F-functions  for  various  grid  densities  are  shown  in  Fig.  12.  The  F- functions  calculated  using  Whitham’s 
method  continue  to  distort  as  the  grid  is  refined  because  the  denominator  in  Eq,  (18)  is  effected  by  Ax,  as 
was  shown  in  Eq.  (23).  Zoomed-in  F-functions  near  the  front  of  the  aircraft  which  includes  the  canard  are 
shovr n  in  Fig,  12(c)  and  12(d).  NVhi thurn  s  method  shows  small  oscillations  or  “scallops.**  These  features 
are  caused  by  the  grid  used  to  calculate  the  canard  s  area  and  lift  not  being  aligned  with  the  canard’s  leading 
and  trailing  edges  as  shown  in  Fig.  13.  This  misalignment  causes  the  area  profile  to  be  stepped;  creating 
a  series  of  positive  to  negative  second  derivatives.  Although  the  modified  Lighthil!  method  uses  the  same 
stepped  area  profile,  the  oscillations  in  the  F-function  are  less  severe  because  of  the  lighter  weight  given  to 
these  points  close  to  the  location  of  interest.  Also,  Whitham’s  method  predicts  a  large  F-function  from  the 
engine  inlet,  as  shown  in  Fig.  12(e).  The  size  of  this  feature  increases  with  increasing  number  of  points. 
In  contrast,  modified  LighthiM’s  method  is  much  less  sensitive  to  these  perturbations.  Figure  14  shows  the 
ground  pressure  profiles  associated  with  the  F-functions  presented  in  Fig.  1 2.  WhiLham’s  method  reasonably 
matches  the  modified  Lighthill’s  method  despite  the  oscillations  associated  with  the  lifting  surfaces  and  the 
engines.  1  his  is  because  these  spikes  are  like  delta  functions  and  contain  very  little  acoustical  energy  and 
quickly  dissipate  with  distance. 
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The  results  of  these  tests  are  summarized  in  Pig.  15.  Figure  15(a)  shows  the  maximum  F-function  value 
for  each  grid  size.  Whitham’s  method  did  not  show  F-funetion  convergence  and  validated  the  analysis  of  the 
previous  section  that  concluded  modified  Lighthiil’s  method  is  the  preferred  method  to  be  used  for  aircraft  F- 
function  calculations  in  RapidF.  A  dashed  vertical  line  at  320  points  is  drawn  in  the  graphs  of  Fig.  15  which  is 
the  lowest  number  of  points  used  in  this  study.  In  Fig.  15(b),  the  maximum  pressure  value  after  propagation 
to  the  ground  is  shown.  The  modified  Lighlhill's  method  maximum  pressure  remains  relatively  constant 
after  320  grid  points.  Although  Whilham’s  maximum  pressure  was  also  flat  after  320  points,  it  was  about  5% 
lower  than  modified  LighthiU’s  pressure.  Finally,  the  run  times  for  each  grid  density  is  shown  in  Fig.  15(c) 
on  a  PC  computer  with  a  quad  core  AMD  640 II  processor  and  12  GB  of  memory.  Both  F-functions  were 
calculated  simultaneously  during  a  single  run;  thus  both  methods  have  the  same  run  times.  The  selection 
of  320  points  allowed:  run  times  to  be  less  than  one  minute;  the  maximum  F-function  to  be  with  3.2%  of 
the  value  using  1 ,600  points;  and  produced  essentially  the  same  ground  pressure  profile  as  predicted  using 
higher  grid  density. 
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(a)  Whithanfs  F-function, 


(b)  Modified  LighlhiN's  F-function. 
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(e)  Whit  ham's  F- function  near  nose. 


(d)  Modified  Lighthilfs  F-function  near  nose. 


(e)  Whitham's  F-funefion  at  engine  inlet. 


(f)  Modified  Lighthiirs  F-function  at  engine  inlet. 


Fig*  12  F-functions  for  various  grid  densities. 
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Pressure,  psf 


Fig,  13  Grid  used  to  obtain  canard  lift  and  volume  distribution. 


(a)  Whitham  s  method.  (b)  Modified  Light  hill  method. 

Fig,  14  Ground  pressure  profiles  lor  various  grid  densities. 
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Max.  F— function,  fl 


(  a)  F-fu  notion.  (b)  Ground  pressures. 


(c)  Run  times. 

Fig*  15  Performance  metrics  for  various  grid  densities. 
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TIL  Rapid F  Sonic  Boom  Pressure  Profiles 


In  this  section,  sonic  booms  from  RapidF  are  compared  to  historical  results  including  legacy  software 
based  on  modified  linear  theory,  wind  tunnel  experiments,  and  (light  test  results.  Since  RapidF  is  based  on 
a  first  order  theory,  it  was  known  prior  to  this  research  (hat  it  is  not  the  most  accurate  software  available. 
The  purpose  of  this  section  is  to  demonstrate  that  RapidF  is  sufficiently  accurate,  defined  as  within  10%  of 
experimental  results.  More  specifically  RapidF’s  sonic  booms  should  have; 

1.  The  proper  shape;  RapidF’s  sonic  boom  should  have  the  correct  number  of  shocks  and  the  shocks 
should  be  within  1 0%  of  the  experimental  location 

2.  Minimum  and  maximum  pressures  within  10%  of  the  experimental  results 

3.  Signal  length  should  be  within  10%  of  the  experimental  signal  length 

This  section  shows  that,  although  it  is  low  fidelity,  based  on  first  order  theory,  RapidF  is  able  to  accurately 
replicate  sonic  booms.  Although  RapidF  occasionally  does  not  meet  the  above  criteria;  overall,  it  consistently 
creates  the  correct  shaped  sonic  boom  and  can  be  trusted  to  predict  the  sonic  boom  of  conceptual  designs. 


A.  Comparison  of  RapidF  to  Wind  lYmncl  Results 

For  an  initial  demonstration  ol  RapidF  s  capabilities,  pressure  profiles  were  compared  to  wind  tunnel 
experiments.  Unlike  sonic  booms  from  (light  tests,  wind  tunnel  studies  lake  place  in  a  more  controlled 
environment,  mitigating  the  effects  of  atmospheric  variables  such  as  humility,  winds,  and  turbulence;  thus, 
making  results  more  repeatable.  NASA  TM  X-1236  [10],  authored  by  Carlson,  McLean  and  Shrout  in  1966, 
was  selected  for  this  purpose.  Despite  its  age,  the  report  has  many  desirable  traits.  First  of  all,  it  presents  a 
complete  geometric  description  of  the  wind  tunnel  model,  allowing  implementation  in  RapidF  with  minimal 
assumptions.  Second,  it  presents  results  for  a  variety  of  test  conditions: 

•  Configuration:  with  and  without  nacelles,  and  two  fuselage  shapes 

a  Speed:  M  =1.41  and  2.0 

o  Lift:  CL  =  0.0,  0.1,  and  0.2 

•  Propagation  angle:  ^  =  0, 45,  90  and  180° 

•  Propagation  distance:  2.5, 5.0  and  10.0  body  lengths  away 
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Furthermore,  it  provided  both  experimental  and  theoretical  results,  allowing  verification  by  comparing 
Rapid F  s  output  to  the  theoretical  data  and  validation  by  comparing  to  the  experimental  results, 

A  diawing  of  the  model  used  in  the  NASA  1M  X-1236  wind  tunnel  test  can  be  seen  in  Fie  16,  with 
dimension  in  centimeters.  The  model’s  engines  were  solid  bodies  of  revolution  with  radii  determined  by 
subtracting  the  inlet  capture  area  from  the  engine  nacelle  geometry.  RapidF’s  model  is  depicted  in  Fig.  17. 
There  are  some  minor  differences  compared  to  the  model  used  in  the  wind  tunnel  tests.  RapidF  used  one 
engine  on  the  centerline  with  the  same  area  profile  of  all  lour  of  the  wind  tunnel  model’s  engines.  This  should 
not  cause  any  errors  for  propagations  directly  beneath  or  above  the  aircraft,  but  might  cause  some  disparities 
for  side  propagations  as  shown  in  Fig.  18.  For  simplicity  only  two  engines  are  shown  in  this  figure.  When  the 
aircraft  with  horizontally  displaced  engines  is  sliced  at  the  Mach  angle,  each  engine  can  create  an  its  own  area 
distribution.  Thus,  the  simplified  engine  on  RapidF’s  model  might  cause  some  errors  for  side  propagation. 
Finally,  the  fillet  where  the  wing  trailing  edge  meets  the  fuselage  is  not  modeled  in  RapidF. 

Like  RapidF,  the  NASA  TM  X-1236’s  theoretical  pressure  distribution  was  calculated  using  modified 
linear  theory,  Thu  area  profile  was  calculated  by  slicing  the  aircraft  at  the  Mach  angle,  corrected  for  angle 
of  attack  by  adapting  a  code  designed  for  wave  drag  calculations  [  14j.  This  code  divided  components  into  a 
maximum  of  30  slices,  then  estimated  the  area  at  each  location  to  include  the  boundary  layer  displacement 
thickness.  The  lift  distribution  was  calculated  to  include  the  nacelle  to  wing  interference  effects.  The  lift 


Fig.  16  Drawing  of  wind  tunnel  model  (reproduced  from  NASA  TM  X-1236  [10]). 
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(a)  Planform.  (b)  side  view. 

Fig.  17  RapidF’s  representation  of  wind  tunnel  model  used  in  NASA  TM  X-1236  [10]. 


Single  engine  on  centerline 


Two  engines  horizontally 
displaced 


Fig,  18  Area  distributions  of  one  versus  two  engines  when  propagating  to  the  side. 

profile  was  converted  Lo  an  equivalent  area  via  Kq.  (20).  Then,  using  (he  total  area  distribution,  the  F-funclion 
was  calculated  using  Whitham’s  equation,  Fq.  (18).  The  pressure  profile  was  found  from  the  F-function  as 
specified  in  NASA  TN  D-3082  [15],  the  same  process  as  RapidF.  In  the  following  discussion  and  figures,  the 
solutions  from  this  process  is  referred  to  as  the  “theoretical’'  results. 
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Id  Fig*  19(a),  the  results  are  shown  for  the  model  with  CL  -  0  at  10  body  lengths.  The  pressure  and 
distance  are  non-dimensionalized  as: 


In  general  RapidF  agrees  very  well  with  both  NASA  TM  X-1236's  experimental  and  theoretical  results.  In 
lact  RapidF  does  a  better  job  capturing  the  small  disturbance  at  x  —  0*1  and  the  wing  s  overpressure  from 
x  —  0.2  to  0.4.  Both  NASA's  theoretical  approximation  and  RapidF  have  some  errors  associated  with  the 
nacelles  from  x  —  0.4  to  0.6*  The  authors  of  NASA  TM  X-1236  were  concerned  about  this  disagreement 
and  noted  that  a  correction  needed  to  be  added  to  account  for  the  soldering  around  the  engines.  However* 
this  adjustment  did  not  fully  account  for  the  difference. 

Pi  ess ure  profiles  are  presented  for  (_x  —  0.1  at  10  body  lengths,  with  and  without  engine  nacelles, 
in  I  igs.  19(b)  and  19(c),  respectfully*  In  both  graphs,  RapidF  properly  captures  the  wing's  overpressure* 
demonstrating  that  RapidF  properly  accounts  tor  lift.  The  inlet  shocks  in  Fig.  19(b)  are  good  indicators  on 
how  well  the  software  predicts  the  effects  from  smaller  components*  The  difference  between  RapidF  and  the 
experimental  results  is  minor  compared  to  Carlson  s  theoretical  pressure  profile.  Carlson  performed  the  tests 
without  nacelles  because  he  speculated  that  the  nacelles  were  causing  problems.  However*  as  can  be  seen  in 
big.  19(c),  even  without  the  engine  nacelles,  Carlson's  theoretical  results  poorly  predict  the  pressure  profile 
at  the  rear  of  the  signature*  RapidF  does  a  better  job,  despite  the  exclusion  of  the  wing-fuselage  fillet. 

Results  for  propagation  angles  of  45°  and  180*  at  10  body  lengths  can  be  seen  in  Figs*  19(d)  and  19(e). 
For  45*  the  agreement  is  very  good  despite  the  fact  that  RapidF  models  the  engines  as  one  nacelle  on  the 
centerline.  Examining  the  pressures  for  180°  in  Fig.  19(e),  RapidF  under- predicts  the  wing's  contribution 
from  x  =  0.2  to  0,3  for  CL 0T.  The  cause  of  this  disagreement  is  not  known  and  is  unexpected  because  for 
all  the  cases  thus  far*  RapidF  has  done  an  equal  or  better  job  matching  the  experimental  results  than  Carlson's 
theoretical  predictions*  Therefore,  an  error  in  RapidF's  lift  or  area  calculations  is  ruled  out.  Furthermore,  the 
results  for  45*  show  proper  application  of  the  trigonometry  for  cases  where  cj>  ^  0.  The  authors  speculate, 
that  the  NASA  report  mis-reported  the  coefficient  of  lift  for  this  case  because  as  shown  in  Fig.  19(e),  if 
RapidF's  CL  is  set  to  0.04,  then  all  the  data  is  in  agreement.  The  final  lest  case  is  for  Mach  2.0,  at  9  body 
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te rights  can  he  seen  in  F  ig.  19(f),  which  assesses  how  well  RapidF  accounts  for  Mach  number  changes.  Like 
the  Mach  !  .41  case  (Fig.  19(b)),  RapidF  more  accurately  predicts  the  shock  locations  from  the  nacelles  than 
the  NASA  theoretical  solution.  Also,  RapidF  more  accurately  predicts  the  signal  duration. 

A  summary  of  the  metrics  used  to  assess  to  accuracy  of  RapidF  are  summarized  in  Table  2.  Items  that 
exceed  the  10%  accuracy  goal  are  in  bold.  While  better  agreement  was  desired,  these  results  should  be  put 
in  perspective.  The  wind  tunnel  model  was  small,  10.54  cm  in  length  and  the  engine  nacelles  were  0.142 
cm  in  diameter.  In  fact,  NASA  TM  X-1236,  was  inspired  by  previous  tests  [16]  which  used  a  2.5  cm  long 
model.  NASA  TX  X-1236’s  test  were  scaled  up  from  these  tests  in  an  attempt  to  mitigate  the  issues  with 
manufacturing,  measuring  forces,  vibrations  and  quantifying  the  boundary  layer  effects  of  very  small  models. 
Although  the  10.45  cm  model  was  better,  some  of  these  issues  still  exist.  For  example,  lift  was  estimated 
(not  measured)  by  setting  the  angle  of  attack,  because  the  model  was  too  small  for  a  force  balance.  The 
angle  of  attack  was  set  based  on  the  theoretical  lift  curve  slope  including  an  estimation  of  model  deformation 
under  test  conditions  with  10  psi  of  dynamic  pressure.  The  authors  of  NASA  TX  X-1236  concluded,  “In 
view  of  the  difficulty  of  consiruciing  accurate  small-scale  models  and  the  difficulty  of  evaluating  boundary- 
layer  effects  for  the  models,  there  remains  some  doubt  whether  extrapolation  of  wind-tunnel  results  (based 
on  theory)  offers  any  increase  in  accuracy  of  (light  overpressure  estimates  compared  with  estimates  based 
on  theory  alone  [10]."  Thus,  the  pressure  transducer  data  should  not  be  interpreted  as  the  absolutely  correct 
solution.  Typical  maximum  pressure  recorded  in  a  signature  was  0.035  pounds  per  square  inch.  Since  only 
35  data  points  were  measured  along  the  sonic  boom,  resolution  is  limited  and  the  shock  locations  had  to  be 
interpolated  to  obtain  the  data  presented  in  Table  2.  Finally,  NASA  TM  X-1236  only  tested  at  a  maximum 
ol  10  body  lengths  away.  Modified  linear  theory  is  better  suited  for  far  field  predictions.  Also,  although 
RapidF  used  550  points  along  the  fuselage,  the  fuselage  description  in  the  NASA  report  only  had  21  points. 

I  here  tore,  RapidF  had  to  interpolate  between  the  data  provided,  and  minimum  improvements  in  results 
should  be  expected  from  the  increased  grid  density.  The  conclusion  from  this  case  is  that  RapidF,  using 
modified  linear  theory,  is  capable  using  a  geometric  description  to  predict  the  pressure  profile  of  a  wind 
tunnel  model  with  various  configurations,  lift,  Mach  number  and  propagation  angles. 
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Table  2  RapidF  percent  difference  from  experimental  data  [10] 


Case 

Shock  Loc.  I 

f  Shock  Loc,  2 

Shock  Loc.  3 

Pm.ni 

Pmin 

Duration 

CL  =  0 

3.8 

-4.0 

N/A 

19 

16 

-5.6 

Cl  =  OJ 

3.5 

-2.3 

0.5 

-1.4 

27 

-1.6 

No  engines 

1.1 

3.7 

N/A 

0.4 

-30 

-03 

<t>  =  45" 

-6.1 

5.3 

N/A 

7.4 

-3.3 

-0.4 

4>  =  180" 

7.8 

N/A 

N/A 

27 

-L5 

0.6 

M  =  2.0 

15 

-3.5 

-0,8 

-10 

-6.7 

-2.0 
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(a)  M  =  1.41,  CL  =  0,0  =  0. 


(c)  M  =  1.41,  Ci,  =  0.1,0  =  0,  without  nacelles. 


<e)  M  -  1.41,  Cl  =0.1,0=  180". 


(b)  M  =  1.41,  CL  =  0.1,0  =  0. 


<d)  M  =  1.41,  Cl  =0.1.0  =  45°. 


(0  M  =  2.0,  CL  =  0.1, 0  =  0. 


Fig.  19  Pressure  profiles  of  NASA  TM  X-1236  [10], 
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B.  Flight  Test  Comparisons 

/.  Lockheed  F- 104  Star  fighter 

Having  demonstrated  the  ability  of  RapidF  to  predict  sonic  booms  in  the  controlled  wind  tunnel 
environment,  a  demonstration  ol  its  ability  to  replicate  sonic  booms  from  flight  tests  is  now  appropriate.  In 
addition  to  the  more  complex  flow  field  with  winds,  turbulence  and  humidity;  replicating  sonic  booms  from 
flight  test  is  also  more  challenging  because  the  test  aircraft  is  much  more  complex  than  a  wind  tunnel  model. 
An  aircraft  has  many  minor  protuberances  such  as  lights,  antennas,  sensors,  trim  tabs,  etc..  Thus,  in  order 
to  properly  predict  its  sonic  boom,  RapidF  must  account  for  the  dominate  features  to  obtain  correct  area 
distribution  and  F-function. 

An  example  of  the  accuracy  of  RapidFs  cross-sectional  area  estimate  is  can  be  seen  in  Fig.  20,  where 
the  cross  sectional  area  of  an  F-104  from  RapidF  is  compared  to  published  data  f  1 7],  It  is  important  to  note 
that  this  is  the  actual  cross-sectional  area,  sliced  perpendicular  to  the  longitudinal  axis,  not  at  the  Mach  angle. 
Although  there  is  slight  differences,  RapidF’s  area  profile  lias  the  correct  shape  and  magnitude.  Thus,  RapidF 
has  enough  fidelity  to  capture  the  dominant  features  of  an  aircraft. 

Figure  2!  compares  RapidF’s  F-function  for  a  15,000  lb  F-104  dying  at  20,000  ft  at  Mach  1.4  to  that 
published  by  Hayes  in  1968  [18].  In  order  to  get  the  lobes  to  line  up,  Hayes’s  F-function  had  to  be  shifted  5 
ft  to  the  left.  By  doing  this  adjustment,  one  can  see  that  both  F- functions  have  the  same  general  shape  for  the 


Fig.  2(1  F-104  cross  sectional  area  distributions  from  RapidF  compared  to  NASA  TR  R-198  [17]. 
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Fig.  21  F-functions  from  RapidF  and  NASA  SP-180  for  a  IS, 000  lb  F-104  cruising  a!  Mach  1.4  and  20,000  ft  [18J. 

first  40  ft,  but  she  magnitudes  differ.  In  order  to  investigate  these  discrepancies,  the  F-function  as  presented 
in  Eq.  (4)  can  be  inverted  to  find  the  area  distribution  [  19|: 

s«i{y)  =  4  /  F(QV(y  -  0  dc  04) 

Jo 

The  area  determined  from  the  above  equation  for  Hayes’s  F-funclion  is  compared  to  that  used  by  RapidF  in 
22.  1  hese  results  differ  from  the  cross  sectional  area  distribution  shown  in  Fig.  20  in  that  they  include 
lift,  boundary  layer  growtii  and  engine  effects.  They  also  differ  in  that  the  aircraft  is  cut  at  the  Mach  angle 
corrected  for  angle  of  attack.  Both  area  profiles  have  the  same  shape  except  that  Hayes’s  area  is  spread  over 
a  larger  distance.  The  effective  length  of  the  aircraft  was  calculated  to  determine  which  area  distribution  is 
correct.  The  component  that  is  farthest  aft  on  the  F-104  is  the  horizontal  tail,  which  terminates  at  54,8  ft  120] 
as  shown  in  Fig.  23.  Since  the  horizontal  tail  is  7.2  ft  above  the  centerline,  its  influence  will  affect  the  area 
distribution  until  61.9  ft.  Thus,  the  effective  length  of  the  F-104  is  61.9  ft  and  after  this  location  the  area 
should  remain  constant.  A  vertical  line  is  drawn  at  this  location  in  Fig.  22.  Since  Hayes’s  area  distribution 
continues  to  change  past  this  location,  it  is  concluded  that  there  must  be  an  error  in  his  area  profile.  As 
another  check,  using  the  same  process,  the  canopy’s  influence  should  be  a  maximum  at  13.7  It,  yet  Hayes’s 
area  distribution  shows  it  at  20  ft.  Therefore,  it  is  concluded  that  our  F-104  F-function  is  more  accurate  than 
Hayes’s  F-function. 
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Fig.  22  F-104  equivalent  area  distributions  from  RapidF  and  NASA  SP-I80  for  same  ease  as  Fig.  21  [18]. 

As  a  check  of  (he  entire  RapidF  system.  Fig.  24(b)  shows  the  sonic  boom  for  an  F-104  Hying  at  Mach 
1.5  and  28,000  ft  compared  to  flight  lest  data  from  NASA  TN  D-2539  [21],  Since  only  surface  weather 
was  published,  RapidF  used  standard  day  conditions  and  no  winds.  Despite  the  limited  information,  RapidF 
properly  predicted  that  the  sonic  boom  is  a  multi-shock,  saw-toothed  signature.  The  shocks  for  the  (light 
test  data  are  considerably  more  smoothed  out  and  RapidF  over-predicts  the  peak  pressures  by  14%  and  13%, 
which  can  be  partially  explained  by  RapidF  under-predicting  the  effects  of  molecular  relaxation.  Addition¬ 
ally,  like  the  wind-tunnel  test  presented  earlier,  the  reader  is  reminded  that  the  experimental  data  has  error 
and  is  not  perfect.  Although  the  authors  pf  NASA  FN  D-2539  do  not  specify  the  uncertainty  of  their  pressure 
transducers,  they  do  state  that  the  have  a  frequency  response  of  only  5,000  Hz,  which  will  cause  shocks  to 
be  more  rounded  and  have  a  lower  peak  pressure  than  reality.  Thus,  the  differences  in  shock  shape  are  of 
little  concern.  Rased  on  this  information,  it  can  be  concluded  that  the  interior  shocks  are  caused  by  the  inlet 
and  wings,  respectfully.  In  summary,  RapidF  estimated  the  interior  shock  locations  within  0.55%  and  10%, 
respectfully  and  predicted  the  duration  within  0.52%  of  the  (light  test  results. 
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Fig.  23  Side  view  of  F-104  showing  influence  of  canopy  and  horizon  la  l  tail  118,20]* 
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(a)  Profile.  (b)  Sonic  boom. 

Fig.  24  Sonic  boom  from  Rapid  F  of  a  F- 104  flying  at  Mach  1,5  and  28,000  ft  compared  to  bight  test  results  [21 J* 


2.  Shaped  Sonic  Boom  Demonstrator 

A  second  example  of  RapidFs  sonic  boom  prediction  capability  is  its  ability  to  replicate  the  sonic  boom 
pressure  profile  from  the  Shaped  Sonic  Boom  Demonstrator  (SSBD),  shown  in  Fig.  25  [22].  Figure  26  shows 
results  for  the  SSBD  flying  a l  Mach  1 .4  and  32,700  ft  [23,  24].  For  this  case  atmospheric  data  was  available 
and  entered  into  RapidF.  The  CFD  solution  presented  was  accomplished  prior  to  (light  by  Lockheed-Martin 
using  S  PLITFLOW  3D  using  an  in  viscid  solver  with  5.5  million  cells.  The  CFD  solution  was  propagated  to 
the  ground  using  PCBoom  [25].  CFD  predicts  a  weak  shock  at  0.02  seconds  that  does  not  exist  in  either  of 
RapidF’s  or  High!  test's  data.  As  summarized  in  Table  3,  RapidF  is  nearly  as  accurate  as  the  CFD-with  all 
metrics  within  the  10%  target,  except  the  interior  shock  location.  CFD  under  predicts  this  location  by  55% 
and  RapidF  over  predicts  it  by  31%.  Overall,  RapidF  properly  predicts  all  other  metrics  within  10%  of  the 
flight  test  results.  Considering  the  speed  of  RapidF,  20  seconds  on  a  PC  for  this  case,  and  the  fact  that  it 
requires  no  mesh  generation,  the  slight  decrease  in  accuracy  of  RapidF  is  an  excellent  trade-off  for  the  time 
savings. 
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Fig-  26  SSBD  sonic  booms  For  Mach  1,41  and  32,700  ft  [23, 24], 

In  previous  publications  [1,  2],  RapidF  was  unable  to  match  the  duration  of  the  SSBD  sonic  boom. 
Fig.  27  shows  the  pressure  profiles  From  Refs,  11,2].  At  the  time  of  publishing  these  papers,  RapidF  used 
LighthilFs  F- function,  Eq.  (28),  and  not  the  modified  LighthilFs  equation,  Hq.  (33).  Furthermore,  RapidF 
converted  the  total  equivalent  area  to  a  radius:  72(C)  —  y/ Sc.q(Qfir.  As  discussed  in  the  previous  sec¬ 
tion,  Ligh th ill's  weighting  double  corrects  area  changes  back  to  centerline,  causing  shorter  than  desired 
F- functions.  Likewise,  the  radius  based  on  total  equivalent  area  also  shortens  the  F-function,  Another  im- 

Table  3  RapidF  and  CFO  percent  error  compared  to  SSBD  Might  test  [23, 24] 


IS  PR  Shock  Loc.  pmax 

Pmin 

Duration 

RapidF 

-53 

31  -10 

"4.1 

-3,7 

CFD 

0.1 1 

-SS  3.9 

8.4 

0.28 
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Fig*  27  Previously  published  [ 1 ,  2]  RapidF  SSBI>  sonic  boom  prediction  [23,  24). 
provement  to  CUBoom  that  cart  be  seen  in  Fig.  26  is  that  the  program  now  estimates  shock  rise  times* 

3.  North  American  XB  70  Valkyrie 

Newly  refined  XB-70  sonic  boom  data  was  published  by  Maglieri,  Henderson  and  Tinelti  in  201  I  [26]. 
Although  the  testing  was  conducting  in  1966  the  data  was  never  fully  reduced  to  a  usable  form.  Maglieri 
recognized  that  the  data  was  valuable  lor  code  verification  because  the  sonic  booms  were  recorded  in  the 
upper  atmosphere,  free  from  the  majority  of  the  distortions  created  by  the  earth’s  boundary  layer;  thus  the 
data  was  “clean.'*  Also,  the  sonic  booms  had  multiple  shocks  and  were  not  simple  N- waves.  Furthermore, 
there  were  multiple  data  runs  at  different  ranges  and  angles,  thus  allowing  many  comparisons*  For  these  tests 
the  XB-70's  weight  ranged  front  324,000  to  350,000  lbs,  the  Mach  number  from  1 .47  to  150  and  altitude 
from  37,584  to  37,917  ft.  The  1966  in-flight  data  was  obtained  by  flying  a  F-104  through  the  sonic  boom 
with  a  specially  designed  pressure  probe*  Due  to  the  dif  ficulty  in  transforming  the  sonic  boom  data  to  a  sLatic 
frame  of  reference,  in  flight  sonic  booms  are  presented  with  respect  to  a  non-dimensional  distance  (distance 
/  signal  length).  The  pressure  transducer  used  had  3%  uncertainty  and  the  reflection  factor  was  estimated  to 
be  accurate  to  within  10%;  yielding  an  overall  uncertainty  in  the  experimental  data  of  10.4%  [26]. 

Images  of  the  XR-70  can  be  seen  in  Fig.  28.  RapidFs  XB-70  model  was  made  from  geometric  descrip¬ 
tion  from  NASA  CR-2GI 1-217078  [27J  and  NASA  TP  1516  [281  and  can  be  seen  in  in  Fig.  29*  The  XB-70 
had  several  features  that  made  inputting  into  RapidF  difficult.  First,  the  XB-70  had  a  hinged  wing  allowing 
the  outer  portion  to  deflect  down  during  flight,  as  shown  in  Fig.  28(c).  This  was  accounted  for  by  projecting 
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the  wing  tip  onto  the  wing  horizontal  plane,  which  appears  as  a  wing  sweep  change  in  Fig,  29(a),  The  ver 
deal  projection  was  modeled  as  ventral  fins,  as  seen  in  Fig,  29(b),  The  canard  and  elevator  deflections  were 
published  for  each  test  ease.  The  canard  data  was  inputted  into  RapidF,  hut  the  elevator  deflections  were  not 
modeled.  Instead  of  using  a  simply  description  of  the  fuselage  as  shown  in  Table  1 ,  a  radius  versus  distance 
data  file  was  created  from  XB-70  drawings,  which  allowed  a  higher  fidelity  mode!  of  the  fuselage.  However, 
the  XB-70  had  a  hinged  nose  similar  to  that  of  the  Concorde.  The  NASA  flight  test  data  was  taken  with  the 
nose  in  the  down  position;  however,  RapidFs  fuselage  model  was  for  the  nose  up  case.  Although,  the  XB-70 
has  two  inlets  which  distribute  airflow  to  six  engines,  the  engines  were  modeled  as  two  large  engines.  The 
inlet  capture  area  is  published,  as  well  as,  the  percent  of  spillage.  This  allowed  the  authors  to  quantify  the 
effective  area  of  the  engine  inlets.  RapidF  models  engines  as  cylinders,  but  the  XB-70’s  engine  pod  has  a 
significant  area  increase  as  shown  in  Fig.  28(a).  This  was  accounted  for  by  adding  area  to  the  fuselage  at  the 
correct  longitudinal  distance  (area  increase  can  be  seen  in  Fig,  29(a)  from  120  to  170  ft).  Furthermore,  the 
inlet  ducts  have  a  system  of  bypasses  to  bleed  off  excessive  airflow  through  the  inlet,  shown  in  Fig.  28(a). 
Although  this  bleed  does  contribute  the  effective  cross  sectional  area,  it  was  not  modeled. 

Sonic  booms  for  26  body  lengths  (BLs)  are  shown  Fig.  30.  Data  was  obtained  by  the  F-104  as  it 
flew  4,727  feel  below  the  XB-70  with  minimal  offset  (<f>  =  0.4^).  RapidF  predicts  two  shocks  from  the 
nose  and  canards  separated  by  30  milliseconds,  matching  the  experiential  results  quite  well,  predicting  the 
overpressure  to  within  5%.  Also,  by  examining,  Fig.  30(a),  one  can  conclude  that  the  two  leading  shocks 
are  caused  my  the  nose  and  the  canards,  respectfully.  Looking  further  down  the  pressure  profile,  RapidF’s 
solution  is  missing  the  interior  shock  at  approximately  0.1  seconds.  This  missing  shock  was  determined  to  be 
caused  by  RapidF  under-predicting  the  effects  of  the  engine  inlets,  as  shown  in  Fig.  30(a).  This  was  rectified 
by  increasing  the  percent  spillage  from  the  published  38%  to  99%.  The  resuks  from  this  modification  is 
shown  on  Fig.  31.  This  example  demonstrates  the  power  of  RapidF.  By  examining  the  F-function  overlaid  on 
the  aircraft  (Fig,  30(a))  and  correlating  portions  of  the  pressure  profile  with  aircraft  components,  the  cause 
of  the  discrepancy  was  quite  intuitive. 

Figure  32  shows  sonic  boom  pressure  profiles  for  an  off-axis  case  with  a  propagation  angle  is  57°.  The 
pressures  were  obtained  2,850  feet  below  and  2,900  feet  to  the  side  of  the  XB-70.  Agian,  RapidF  used  99% 
spillage  This  example  is  more  technically  challenging  than  the  previous  example  because  it  is  closer  to  the 
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(a)  Plunform, 


(b)  Side  view. 


Projection  onto 
wing’s  plane 


Wing  tip  in 
65°  position 


Projection 
onto  vertical 
plane 


(c)  From  view  wilh  wintip  projections. 


Fig.  2K  XB-70  drawings,  reproduced  from  RoL  [26]* 


(a)  Plan  form. 


<b>  Side  view. 


Fig,  29  RapidF  XB-70  model* 


aircraft,  19  BLs;  thus,  the  pressure  protile  is  more  complex,  RapidF  properly  predicts  the  shock  locations 
and  magnitudes,  except  for  the  trailing  shock  which  is  over- predicted  by  43%. 

An  over-thc-top  case  is  presented  in  Fig.  33.  In  this  case  the  F- 1 04  flew  3,290  feet  above  and  7,1 00  feet 
to  the  side  of  the  XB-70,  for  a  propagation  angle  of  114°.  In  Fig,  33(a),  the  pressure  profile  is  shown  for 
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F-funct.  x  200  ft1  2  F-funct.  x  200  ft 


(a)  PI  an  form  and  F-fu  notion. 


Distance/Signal  Length 


(b)  Sonic  boom. 


Fig,  30  Sonic  booms  for  XB-70  at  20  body  lengths  [26], 
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(a)  Plan  form  and  F-  function. 


(b)  Sonic  boom. 


Fig,  31  Sonic  booms  for  XB-70  at  26  body  lengths  with  99%  Spillage  [26], 


Distance/Signal  Length 


Fig,  32  XB-70  sonic  boom  for  <j>  =  57  and  19  BLs  [26], 
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(a)  Sonic  Boom, 


tb)  No  engines  sonic  boom. 


Fig,  33  Over-the-top  sonic  boom  For  XB-70  at  42  body  lengths,  <p  =  115°  [26], 


(a)  Example  of  slicing  from  above,  (b)  Engines  placed  above  wings. 

Fig,  34  Wing  shielding  engines, 

the  standard  RapidF  XB-70  model.  The  interior  shock  caused  by  the  wings  and  inlets  is  over-predicted  by 
173%.  This  is  because  when  RapidF  slices  the  aircraft  along  Mach  lines,  it  includes  the  fuselage,  wings  and 
engines,  as  shown  in  Fig,  34(a).  Thus,  the  engine  inlets  should  contribute  significantly  to  the  sonic  boom. 
However,  the  wing  serves  as  a  barrier  or  shield.  This  prevents  the  engine  from  contributing  to  the  sonic 
boom  when  the  engine  is  in  the  shadow  of  the  wing.  This  concept  has  been  used  in  reverse  as  a  sonic  boom 
minimization  strategy.  In  order  to  reduce  the  sonic  boom  on  the  ground,  the  engines  are  mounted  over  the 
wings  1 29,  30],  1  his  configuration  allows  the  downward  inlet  shocks  to  be  reflected  upward  by  the  wings,  as 
shown  in  big,  34(b).  Thus,  in  order  to  improve  the  accuracy  of  the  RapidF  pressure  profile,  for  over-the-top 
sonic  booms,  the  engine  is  omitted  from  the  input,  with  much  improved  results  as  shown  in  Fig.  33{b), 

Finally,  the  ground  pressure  profiles  are  shown  in  Figs,  .35,  Pressure  data  was  taken  from  arrays  of 
microphones,  Each  array  had  six  microphones.  The  accuracy  of  the  system  was  estimated  to  be  11%  by 
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(a)  Plan  form  and  F-function. 


(b>  Sunk  boom. 


Fig*  35  Sonic  booms  for  XB-70  at  the  ground,  185  body  lengths  [26]* 


comparing  data  from  al!  arrays.  Since  the  ground  measurements  did  not  have  the  issues  with  a  moving 
reference  frame,  the  results  in  Figs.  35  are  presented  with  a  dimensional  horizontal  axis  of  seconds.  RapidF 
properly  predicts  the  shock  locations  and  magnitude.  However,  it  over  predicts  the  duration  by  15%.  For 
this  example,  RapidF  used  the  published  XB-70  inlet  spillage  values.  Since  the  spillage  only  needed  to  be 
modified  for  the  data  taken  from  19  to  26  BLs,  the  authors  conclude  that  RapidF’s  simple  engine  model  is 
sufficiently  accurate  for  far  field  sonic  boom  pressure  profile  predicts. 

The  results  for  all  XR-70  cases  arc  summarized  in  Table  4.  The  most  significant  errors  are  for  the 
minimum  pressure  on  the  off  axis  cases.  Since  RapidF  was  able  to  properly  replicate  the  trailing  shock  for  all 
previous  examples,  to  include  off  axis  cases  such  as  Fig.  19(d);  the  authors  conclude  that  these  large  errors 
are  unique  to  the  XB-70  and  are  probably  due  to  the  abrupt  termination  of  all  components  (wings,  fuselage, 
engines  and  vertical  tails)  at  189  feet,  as  shown  in  Fig.  28.  Although  some  other  metrics  exceed  the  10% 
target,  since  the  accuracy  of  the  experimental  data  is  on  the  order  of  10%,  these  differences  are  not  significant. 
The  overall  success  ol  RapidF  predicting  the  pressure  profiles  for  a  complex  aircraft  like  the  XB-70  indicates 
its  reliability.  Despite  the  fact  that  RapidF  did  not  model  the  XB-70' s  nose,  ailerons,  and  bleed  air  system. 
RapidF  did  model  the  dominate  features  of  the  aircraft  which  enabled  it  to  replicate  the  flight  Lest  data. 

IV.  Conclusions 

The  intent  of  this  paper  was  to  demonstrate  the  accuracy  of  RapidF  to  replicate  historical  sonic  booms, 
so  that  the  user  can  trust  the  software  to  predict  the  sonic  boom  of  his  design.  RapidF  is  successful  at  this  task 
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Table  4  Percent  error  of  Rapid  F  compared  to  XB-70  bight  test  [26] 


Case 

Shock  Loc, ! 

Shock  Loc.  2  Shock  Loc.  3 

Ptrla.  r 

Pt  rt  i  ti 

Duration 

BL  =  26, <f>  =  0.4" ,  Sp  =  99% 

-7.2 

15 

9.6 

4.9 

- 

BL  =  19,  <j>  =  57",  Sp  =  99% 

-17 

0.6  -8.7 

16 

43 

- 

BL  =  189,  <p  =  0",  Sp  =  38% 

4,5 

- 

-3.3 

-1.9 

15 

BL  =  42,  ^  =  115”,  Sp  =  0% 

6.6 

-3,3 

16 

61 

- 

because  it  uses  LighthilFs  F-f unction  equation  which  was  shown  to  be  superior  to  that  of  Whitburn,  because 
it  converges  as  grid  density  increases  and  because  it  can  accommodate  discontinuities  in  area  profiles .  The 
weighting  in  LighthilFs  equation  was  modified  since  RapidF  already  corrects  for  thickness  when  it  slices  the 
aircraft  to  obtain  an  ax i-sym metric  area  profile.  RapidF,  based  on  modified  linear  theory,  was  able  to  predict 
sonic  boom  pressure  profiles  lor  wind  tunnel  an  flight  tests.  However,  the  XR-70's  complicated  geometry, 
unique  nacelle  design  and  the  bleed  air  made  st  difficult  to  predict  its  sonic  boom.  Despite  these  challenges,  it 
was  able  to  predict  the  correct  shape  and  duration  for  multiple  cases.  In  general,  RapidF  and  modified  linear 
theory  arc  capable  of  producing  a  first-order  pressure  profile,  as  long  as  the  aircraft  model  can  be  resolved 
into  an  axi-symmetric  area  profile.  With  run  times  less  than  a  minute,  it  is  ideal  For  preliminary  design. 
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